ALGEBRAIC STRUCTURE OF QUASIRADIAL 
SOLUTIONS TO THE 7-HARMONIC EQUATION 

Vladimir G. Tkachev 

We obtain an explicit representation for quasiradial 7-har- 
monic functions, which shows that these functions have es- 
sentially algebraic nature. In particular, we give a complete 
description of all 7 which admit algebraic quasiradial solu- 
tions. Unlike the cases 7 = 00 and 7 = 1, only finitely many 
algebraic solutions is shown to exist for any fixed | —y | > 1. 
Moreover, there is a special extremal series of 7 which exactly 
corresponds to the well-known ideal m-atomic gas adiabatic 
constant 7 = ^P^. 

' 2m-\-l 

1. Introduction 

We study specific solutions to the following quasilinear equation 
(1) 

u xx ((7 + l)u 2 x + (7 - l)Uy) + Au xy u x u y + u yy ((j + l)Uy + (7 - l)u 2 x ) = 

where I7I > 1 or 7 = 1. Let L^[u] denote the left-hand side of this equation. 
A solution of the form 

(2) U (x,y) = p k f(9), k>l, 

where p and 6 are the polar coordinates in the (x, y)-plane, is said to be a 
quasiradial. The origin of this study goes back to the well-known p-Laplace 
equation 

div(|Vu| p - 2 W) = 0, 

which is the divergence form (1) with p = -^y. We call solutions to (1) 
7-harmonic functions. 

The existence and integral representations for quasiradial 7-harmonic 
functions were previously established by G. Aronsson in [l]-[5]. In particu- 
lar, it was shown in [2] (see also [15] and [3]) that single-valued quasiradial 
7-harmonic functions do exist only for those exponents k in (2) which satisfy 
the characteristic equation 

(3) (27V-l)(7 + l)A; 2 -2(iV 2 7 + 2iV-l)A; + iV 2 (l + 7) = 0, n G N. 

We refer to the corresponding solution as the N-solution to (1). One can 
easily see that the 0-solutions are constants, and the 1-solutions are linear 
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functions. In what follows, such solutions are said to be the trivial solutions 
to(l). 

The limit case 7 = 00 reduces to the standard Laplace equation Au = 0, 
and the corresponding iV-solutions are harmonic polynomials of degree N. 
Note that harmonic polynomials are algebraic functions, i.e. they satisfy 
some (actually, trivial one) polynomial identity P(x, y, u) = 0. We show 
that this property is still valid for iV-solutions of the Aronsson equation (i.e. 
7 = 1) 

(4) U'xxU'x ^Uxy^x^y ~t~ ^yy^y — 0, 

It is the aim of this paper to study this phenomena for general 7's. To 
make this point more explicit, we have to note that for 7 7^ 00 a (weak) solu- 
tion of (1) is normally in the class C 1,a . In particular, quasiradial solutions 
have a Holder singularity near the origin, and one should consider them as 
'singular solutions' (the terminology is borrowed from [4]). This non-regular 
character is a consequence of the general situation for 7-harmonic functions 
near their singular points (i.e. at the points at which |V«| = 0), see [17], 
[8], [11]. 

Our key result is the following explicit parametric representations for Ab- 
solutions: 

x + l y = e*HK\C\ 2{N - 1) + C 2N - 1 ), 
( } u N = C\(\ k ( 2N -^- N -Re( N . 

Here £ £ C is the parametrization variable, (j) is an arbitrary constant, 
k = k(N,^) is the biggest root of (3), and fi is defined by (41) below. We 
show also that (5) represents an entire graph over the (x, y)-plane when ( 
runs the complex plane C. 

In particular, it immediately follows from (3) and (5), that is an 
algebraic function whenever k(N, 7) is a rational number. This makes more 
explicit the mentioned above Holderian behaviour of quasiradial solutions 
at their singular points. 

In contrast with (4), we show that for any rational [7) > 1, the class 
of algebraic TV-solutions is necessary finite in the sense that the following 
upper estimate holds N < 2 p 2 > wnere 7 = p/q, with q and p to be 
co-prime numbers, and \_x\ denotes the integer part of x. In particular, this 
yields the absence of algebraic solutions for all integer numbers 7, [7) > 2. 

Denote by A the set of all 7 € Q, [7! > 1, such that (1) admits nontrivial 
algebraic iV-solutions. Then 7 G A iff —7 <G A (see Section 4 below). In 
Section 5 we show that for all rational 7 = p/q G A, 7 > 1, the following 
bilateral estimate holds 



2 + q<p<q 2 -2. 



ALGEBRAIC STRUCTURE OF QUASIRADIAL SOLUTIONS 



3 



The last inequalities are sharp. Moreover, in Section 6 we prove that 
equality q = p + 2 holds iff q is an odd number, q > 3. This yields the 
so-called minimal series 

(6) ^=^TT " = 2,3,4,..., 

and N is the index of the corresponding (a unique for the given 7) algebraic 
quasiradial solution. 

We end this introduction with one possible physical interpretation of (6), 
which gives a motivation of our choice of 7 instead of p. Namely, observe 
that (1) is a homogeneous form of the gas dynamics equation 

(7) L 7 [</>] = 

for the potential of the gas velocity <fi [7, p. 9]. The parameter 7 in (7) 
is the so-called adiabatic gas constant, that is the ratio of the gas' specific 
heats at constant volume and constant pressure (see, e.g. [14]). Note that, 
for all known gas models due to a specific combinatoric nature of the adia- 
batic constant 7, it can be evaluated in terms of the freedom degrees of the 
corresponding gas. In particular, it follows that 7 a rational number. 

Moreover, the most important for applications is the simple gas which 
consists of m atoms, 7 is given by the following ratio 

2m + 3 

(8 7m =0 TT' m = l,2,..., 

2m + 1 

(e.g., m = 2, or 72 = 7/5, describes the standard Earth's atmosphere). 

Note that eq. (7) is of non-degenerate elliptic type for all adiabatic ex- 
ponents 7 > 1. In the recent paper of I. Zorina and the author [16], it is 
proved that for any integer N > 2 there exists a solution of (7) with non- 
trivial polynomial growth > 1 (where kjy is defined by (3)) which is a 
real analytic function in the whole M. 2 . 

Now, the quasiradial iV-solutions to (1) can be naturally regarded as the 
cones (after a suitable scale renormalization) over the corresponding N- 
solutions to (7). In other words, (1) represents a microscopy level of the gas 
flow. In this connection, the coincidence of the minimal series (6) and the 
natural adiabatic constants (8) implicitly underlines the essence of algebraic 
character of the corresponding ^-solutions. We observe also, that in this 
case we have N = m, i.e. the atomic number is equal to the index of the 
corresponding iV-solution. 



2. The separation equation 

In this section we study the basic properties of the wave function f(9). For 
technical reasons, we assume that I7I > 1. The case 7 = 1 requires a few 
more care because of degenerate character of the separate equation (13). 
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However, all the formulated below results are still valid in this limit case if 
we suppose that k > 1. 

The separation of variables in (2) yields the following ordinary differential 
equation 



(9) / "(( 7 -l)fc 2 / 2 + (7 + l)/' 2 )+ 



+ 



/ /2 (fe(7 + 3)-2) + ((l + 7 )fe-2)A : 2 / 2 



fk = 0, 



where the prime denotes the derivative with respect to 9. Letting W 
f' 2 (9), Z = f 2 (9), we can rewrite (9) as 

dW _ fc((7 + 3)fc - 2)W + fe 3 ((7 + l)k - 2)Z 
~dZ ~ (7 + 1)W + (7 - l)k 2 Z ' 

which splits into the following linear system 

W'(0 = -k{( 1 + 2,)k-2)W -k\( 1 + l)k-2)Z 

[W) z'(0 = ( 1 + i) W + { 1 -i)k 2 z. 

One can easily verify that 

W(0 = dk 2 e- 2k ^ + C 2 k((j + l)k - 2)e 2 ( fc " fc2 )« 
Z(0 = _C ie - 2fc2 «-C 2 ( 7 + l)e 2 ( fe - fe2 )«, 



is the general solution of (10), where C\ and C2 are arbitrary constants. 
Then 

2k 

(11) W + k 2 Z = -kr, k ~ l C 2 W + \ 2 Z = -V h Ci, 

7+1 

where r\ = e~ 2 ^ k and A 2 = k 2 — ^j. The right-hand side of the latter identity 
is positive for all \j\ > 1 and k > 1, so we can define 



(12) X:=Jk 2 2k 



7 + 1 

Then elimination of rj in (11) yields (W + k 2 Z) k = C 3 (W + X 2 Z) k ~ 1 . Since 
(9) is a homogeneous equation, it suffices to study only the case C3 = 1. 
Thus, we have the following first order differential equation 

(13) (f' 2 (9) + k 2 f 2 (9)) k = (f 2 (9) + \ 2 f 2 (9)) k -\ 

We introduce new phase variables z = f(9), w = f'(9), and define the set 

(14) r = {(z, w) G R 2 : (w 2 + k 2 z 2 f = (w 2 + X 2 z 2 ) k ~\ w 2 + z 2 + 0}. 

Observe that the intersection T with the Oz-axis consists of exactly two 
points (the apexes): = (±zo,0), where 

(15) z = X k ~ 1 k~ k . 
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Lemma 2.1. Let \j\ > 1 and k > 1. Then T a real analytic closed Jordan 
curve. Moreover, F\ {A + ,A~} splits in two mutually symmetric graphs 
(which have no common points with Oz). 

Proof. The first statement easily follows from the representation of T in the 
polar coordinates z = r cos a, w = r sin a: 

(sin 2 a + A 2 cos 2 a)^" 1 )/ 2 

v = . 

(sin 2 a + k 2 cos 2 a) k / 2 

Now, consider V as the 0-level set of the function 

(16) F(z; w) = (w 2 + k 2 z 2 ) k - (w 2 + A 2 * 2 )* -1 . 

We claim that / on T \ {A + , A~}. Indeed, suppose F( n (zi,wi) = 
and wi / 0. Then we have 

(17) ^-F' w { Zl , Wl ) = k(w 2 + k 2 z 2 ) k ~ l - (k - l)(w 2 + A 2 z 2 ) fc " 2 = 0, 

which together with F(zi,w\) = and I7I > 1 implies 

1-7 2 



w 2 = J^fk 2 zl < 0. 

It follows that w\ = z\ = 0, which contradicts the definition of T and implies 
our claim. Thus, T \ {A + ,A~} splits into union of two graphs with respect 
to the Oz-axis and the lemma is proved. □ □ 

Now, we construct a special solution f{6) of (13), satisfying the initial 
condition /(0) = zq. With the above notation we have 

dz f{e)de 

Define 

(i8 » 

where £ € T, and the integral is taken clockwise along the arc (A + ,^) of 
r. The last integrand a priori has singular behavior when w vanishes (i.e. 
for £ = A ± ). But, it can be shown that these singularities are removable. 
Indeed, using again the representation of T as the 0-level set of function (16) 
we find 

<te_ F^dw _ k(w 2 + fcV)^ 1 - (fc - l)(u> 2 + AV) fc ~ 2 dw 
( ' ^ ~ ~~F>T ~ ~ k*(w 2 + k 2 z 2 ) k ~ x — {k — l)X 2 {w 2 + \ 2 z 2 ) k ~ 2 '~' 

Now, to show that (18) has no singularity it suffices only to verify that the 
denominator of the right-hand ratio in (19) is non-zero in a neighborhood of 
the apexes. The corresponding values at A ± are equal A 2 ^ -1 ^^ -2 ^ / 0. 
Therefore, the integral in (18) is well defined, and it follows that 6(£) is an 
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analytic function of £ (in the sense that 0(£(r)) is analytic for any analytic 
parametrization £(r)). 

Next, observe that dz/w > within our convention. Thus, applying 
Lemma 2.1, we conclude that 0(£) is a strictly increasing function when £ 
runs T in clockwise direction. Define function by letting /fc(6(£)) = 

where z(£) is the projection of £ onto the Oz-axis. Clearly, fk(0) 
becomes a real analytic periodic function, which is defined now in R. By 
virtue of the symmetry of T, the one-quoter period T satisfies 

<»> ? =- r - 

where B - = (-1,0) E I\ 
Lemma 2.2. 

fc-1 



(21) T = 2vr ^1 

Proof. Define a new variable t by letting it; = —zt. Then, applying (14) we 
obtain the following parameterization of the arc A + B~ 

z (t) =(t 2 + A 2 )( fc - 1 )/ 2 (t 2 + fc 2 )- fc / 2 

w(t) =-t(t 2 + A 2 )( fc - 1 )/ 2 (t 2 + A: 2 )- fc / 2 , 

When t runs between and +oo the corresponding point £(i) = (z(i), 
runs A + B~ in the clockwise direction. Moreover, we have 



d*(t) =t{t z + \ z ) — (t z + k z )-—\k\k-l)-k\ z -t z j dt, 
which yields 

(23) - (-*— - ±=±-\ dt 

[ ' W{t) ~ \t^ + k 2 t 2 + X 2) at - 

Integration of (23) gives 

I — i i 

(24) @(z(t), w(t)) = arctan — arctan — . 

k A A 

So, letting t -> +oo in (24), we obtain by virtue of (20): \ = § (l - 

and (21) follows. □ □ 

From the definition of T we infer /£(6(£)) = where z(£) is the 

projection of £ onto the Oz-axis. In particular, 

Tr? 

(25) /*(*) = ^ 9 = —, neZ. 

Hence, the constructed function fk(0) satisfies (13) with the initial data 

\ k—i 

/ife(0) = 0, f k (0) = z o = - jpr . 
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It is obvious that (by virtue of the autonomic character of (13)) a general 
solution of (13) has the form f(9) = Cfk{0 + a), where C and a are arbitrary- 
constants. 

Corollary 2.3. Let T be defined by (21). For 9 € (-T/4,T/4) we have the 
following parametrization for fk{9) 

/fc=(t 2 + A 2 ) (fc-l)/2 (t 2 + A .2 ) -fc/2 

(26) t k-1 t 

6 = arctan — arctan — , t £ K, 

k A A 

and for other values of 9, f{9) satisfies the symmetry rules: 

/*(<?) = -/*(f-0)> AH?) = /*(<?)■ 

Our further objective is to characterize all values of k > 1 which support 
the 27r-periodic wave functions fk(9), that is + 2ir) = fk(0). In fact, 
the latter condition is equivalent to absence of multivalued branches of the 
corresponding quasiradial solution of (1). The following assertion is a direct 
consequence of Lemma 2.2. 

Proposition 2.4. Let \j\ > 1 and k > 1. Then fk(9 + 2tt) = fk(9) if and 
only if the following equality holds 

where N denotes the set of all positive integers. In this case fk(9) is a 
2ir / N -periodic function. 

The latter statement is not new. It has appeared in [2] for 7 = 1 and 
for general 7 in [3], [15]. However, our approach to this result seems to be 
complementary to these previous works and allows us to arrive at an explicit 
representation for the quasiradial solutions in the next section. 

Proposition 2.5. Let \j\ > 1 and N £ N be given. Then there exists a 
unique k = £(7, N) > 1 such that (9) admits a 2n / N -periodic solution. 

Proof. The trivial case TV = 1 by (27) gives k = 1. Let N > 2. Then (27) is 
equivalent to the following quadratic equation 

(28) (2N - 1)( 7 + l)k 2 - 2(iV 2 7 + 2N - l)k + N 2 {1 + 7) = 
which has two separate roots because it discriminant is strictly positive: 

D = A(N - l) 2 (iV 2 7 2 - 2N + 1) > 4(iV - l) 3 > 0. 
Then one can easily infer from the Viete theorem that 

(29) ^- 1 )( fc2 - 1 ) = - 1 2 ^f -7Tl <0 < 
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where k\ ^ ki are the roots of (28). Inequality (29) implies k\ < 1 < &2, so 
that exactly one root k-i > 1 is consistent with our constraint k > 1. □ □ 

3. Absolutions 

From now on, we adopt a new notation fx for the iV-th wave function /& 
with k = k(-f,N). 

Definition 3.1. Let N € N. The quasiradial solution of the form 

u N (x,y) := Cp k f N (6), 

where C is an arbitrary constant, is said to be a basic TV-solution of (1). 
Similarly, u = C p k /n(6 + a) with an arbitrary aGlis said to be a (general) 
iV-solution. 

Theorem 3.2. Let \ j\ > 1, N G N, and k = k(j,N) be the biggest root of 
(28). Then the basic N -solution has the following representation 

x = h 2N - l ({k + A) cos t + (k - A) cos(2iV - l)r), 

y = h 2N ~ 1 ((k + A) sinr - (k - A) sin(2iV - l)r), 

u N = Ch k <? N - 1 '* cosTVr, 

where A is defined by (27), and r € [0;27r], h > are i/ie variables of 
parametrization. 

Proof. By virtue of (26) we have the following parametrization for /n(6) 

/jv=(t 2 + A 2 ) (fc-l)/2 (t 2 + fc 2 ) -fc/2 

(30) i fc-1 t 
t/ = arctan — arctan — . 

k A A 

Define a new variable r by 

(31) t = \t*a(NT), re(-^,^). 
Then we have from (30) and (27) 

arctan I 



(32) 6 + (N — l)r = arctan ( £ tan(AV 



which by (31) yields 

(33) i = fctaii(0 + (iV-l)T). 

Inserting (31) and (33) into the first identity in (30) we get 



f N =(t 2 + X 2 )^(t 2 + k 2 )-l 



(34) =z (1 + tan 2 (Nt)) 2 (l + tan 2 (# + (N - l)r)) 

=Zq cos Nt 



cos(# + (iV — 1)t) x * 



cos A^T 
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where zq is defined by (15). On the other hand, 

(35) cos(0 + (N — 1)t) = cos6lcos(iV - l)r • [1 - tan(iV - l)r • tan0]. 
Now, we apply to (35) the addition formula 

1 - tanp • tan(/3 - p) = ■ 1 

cos^ p 1 + tan p tan p 

with p = (N — 1)t, /? = arctan (^ tan(A/"r)) . By virtue of (32), we have 
9 = (3 — p, therefore 

i i 

l-tan(iV - l)rtan6> = 



cos 2 (iV - l)r 1 + A tan(iV - l)r tan Nt 



(36) 



A; cos iVr 



cos(7V - l)r fe cos(iV - 1)t cos iVr + A sin(iV - l)r sin iVr 
1 2k cos iVr 



cos(iV - l)r ' (A; + A) cos r + (A; - A) cos(2iV - l)r ' 
Then, applying (35) and (36) to (34) we obtain 

2k cos 6> 

/JV = Z ° C ° S 7VT [(fc + A)cosr + (A ; -A)cos(2iV-l)r. 
Taking into account that uw(x,y) = p k fN and x = pcos9 we find 

2k cos 



ui\[(x,y) = zq cos Nt 
= (2k) k z 



i k 



_(k + X) cos T + (k-X) cos(2N - l)r _ 

fc 



(A; + A) cos T + (k-X) cos(2iV - l)r _ 

Setting h 2N ~ l for the expression in the last brackets we arrive at 

x = h 2N ~ l ((k + A) cos r + (jfe - A) cos(2AT - l)r), 

and Uiv = C N h k( - 2N ~V cos Nt, where Cat = {2k) k z = 2 k X k ~ l . 

Finally, to express y we eliminate the polar coordinates as follows 

y _ A tan A^r - ktan(N - l)r 

x~ ~~ A; + A tan A/t tan(iV - l)r ~ 

(fc + A) sin r - (jfe - A) sin(2N - l)r 
~ (A; + A) cos r + (k - A) cos(2N - l)r ' 

Thus we get (3.2) for all r G (-tt/2N,it/2N). Using the analyticity of 
/at(#) we conclude that (3.2) is valid for all r. The theorem is proved 
completely. □ □ 
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In order to simplify (3.2) we make use a special intermediate parameter 
\x. Namely, in the above notation, put 

(37) " = —\> 

so A = k{n — l)/(/i + 1), and we have from (12) 

V ^ 2/i( 7 + l) 

Then, an easy computation shows that (27) becomes 

u 2 - 1 

(39) iV = — t-- 

V ' 2(^7-1) 

Now we observe that for 7 > 1 we have k > A > 0, so that fx > 1. Similarly, 
7 < — 1 implies /i < — 1. On the other hand, considering (39) as a quadratic 
equation for 11 

(40) F(n) := li 2 - 2 7 iV> + (2JV - 1) = 0, 

wefindF(l) =2JV(l-7) and F(-l) =2iV(l + 7 ). If 7 > 1 then F(l) < 0, 
i.e. exactly one root of (40) 

fi+ = AT 7 + y/N^f -2N + 1 

agrees the above constrain fx > 1. Similarly, for 7 < —1 we have 



/x" = Nj — y/N^f - 2N + 1. 

Define 



(41) fx^fx^N) 



Nj+ yjN 2 j 2 -2N + 1, for7>l; 
iV 7 - y/iVV -2N + 1, for 7 < -1. 
For the further purpose, notice that 

(42) M (- 7 ,iV) = - M (7,JV). 

Now, by using the homogeneity character of (43), one can rewrite it as 
follows 

x = h 2N - 1 (ncosT + cos(2AT - l)r), 

(43) y = h 2N - 1 (fi sin t -sin(2N -1)t), 

UN = Ch (2N ~^ k cos Nr. 

A general TV-solution can be obtained from a certain basic iV-solution by 
a suitable rotation of in the (x, y) plane. Let 

x 1 = x cos ip + y sin ip, y' = — x sin -ip + y cos tjj, 
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be such a rotation. Then (43) implies 

x = h 2N - 1 (ncosT + cos((2iV - l)r + 2Nip)), 

(44) y = h 2N - l (psiiiT - sin((2iV - l)r + 2Ni/>)), 
u = u N ^ = Ch^ N ~^ k cos N(t + V), 

In particular, ujv,o = un- Thus, (44) gives the representation for general 
N-solutions to (1). 

By putting in (44) X = hcosr and Y = fa sin r, we obtain an algebraic 
representation of a basic iV-solution u(x, y) (see also an equivalent complex 
form (5) given in the Introduction) 

x = fiX(X 2 + Y 2 ) N ~ l + Re(X + iY) 2N ~\ 

(45) y = VY(X 2 + Y 2 ) N ^ - Im(X + iY) 2N ~\ 
UN = C(X 2 + Y 2 ) * Re(X + iY) . 



Corollary 3.3. All N-solutions are quasialgebraic functions in the sense 
that u N is an algebraic function, where a = k ^ 2 N-i)-N ' 

To illustrate the last property, we briefly mention the following well-known 
example. Note (see, also [2] and [10] for further examples) that u 2 = x 4 / 3 — 
y 4 / 3 is a basic 2-solution of (4). Then it is easily verified that u = U2(x,y) 
satisfies the following polynomial identity 

27x 4 y 4 u 3 = (x 4 -y 4 -u 3 f. 



4. Conjugate solutions 

We recall that the main equation (1) can be represented as a p-Laplace 
equation for p = It is well known fact (cf. [5], [6]) that in two- 

dimensional case there is the canonical correspondence between p-harmonic 
and p'-harmonic functions for 

- + i = i. 

p p 

More precisely, given a solution u of equation L 7 [u] = we define the 
conjugate function U by 

(46) U x = \Vu\ 2 ^-^u y , U y = -\Vu\ 2 ^V Ux , 17(0) =0. 

Note that U(x, y) is not necessarily a single valued function. But at least 
locally, U is a quasiradial solution of the conjugated equation 



(47) 



L_ 7 [C7] = 0. 
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It turns out that there is a simple relation between the conjugate quasira- 
dial solutions. We define an adjoint (to basic one) iV-solution follows 

x = h™- 1 ^ cost - cos((2A - l)r)), 



(48) 



y = h 2N ' 1 (fismT + sin((2A - l)r)), 



u* N = Ch^ 2N -^ k sin Nr. 
Applying (44), we obtain an equivalent definition: 

U* N = U N _ n /2N- 

We have also u*£ = —un- The functions un and u* N form a conjugate pair, 
analogous to conjugate harmonic functions. More precisely, if 7 = 1 one can 
easily derive from the above representation that un(x, y) = Re(x + iy) N and 
u* N (x, y) = Im(x + iy) N . 

Theorem 4.1. Let be a basic N -solution of (1) and be an adjoint 
N -solution of (47). Then there is a constant c such that cun and form 
the conjugate pair in the sense (46). 

Proof. Without loss of generality, we can assume N > 2. Let un be an 
Absolution of (1) and U be the corresponding conjugate function defined by 
(46), normalized by U(0) = 0. It follows from homogeneity of un that U is 
a homogeneous function as well. Hence, there is a real (3 such that 

(49) U = r p G(9). 

Here G{6) is a priori a multivalued function. We will prove that U = 
with a suitable constant C in (48). 

First, notice that U is a quasiradial solution of (47). Denote by k the 
growth exponent of un- Then k > 1 and the components of the gradient 
Vun are homogeneous functions of order (k — 1). Moreover, 

NuNl 2 =r 2k-2 {k 2 f 2 N{d) + f ,2 m 

In particular, Vu^ ^ for r ^ 0. Applying (46) gives 
U'J ,_ 



vu 



\Vu N \ 



2/( 7 -l) 



—u 



N,y 
1 



N,x. 



7+1 

7 -i 



Gi{9) 
G 2 {9) 



where Gi{6) are certain 27r-periodic functions of 9. Since | — y | > 1 and k > 1, 

"i + 1 

= (k- 1)J-L- + 1 > 1. 
7 - 1 

On the other hand, let k* be the growth exponent of U^. Let /j, and 
H* be the corresponding auxiliary parameters defined by (37) for k and k* 
respectively. Then it follows that from (42) that \i* = —fi, and from (38) 
we have 

(1 + M) 2 



k 



2(1 + 7 )/i 
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and 

= (1 + ^) 2 = (l-^) 2 
2(1 + 7* V* 2( 7 -l)/x' 

where 7* = —7. Hence 

(50) fc(l + 7) + fc*(l - 7) = fc(l + 7) + k*(l + 7*) = 2. 

Thus, by (50) we have 

/3 = ( *_ 1) l±i + 1 = M2±i>_X±i + 1 = 

7 — 1 7 — 1 7 — 1 

_2-fc*(- 7 + l) 2 _ r 



7 — 1 7 — 1 

The latter means that U is an iV-solution of (47) with (3 = k* (in particular, 
the function G in (49) is a single valued function). 

It remains only to show that U is an adjoint iV-solution. Since U is an 
TV-solution, there is ip such that U = U^^ in representation (44). Thus, we 
have to prove that Nip = ±7r/2 mod it. 

Choosing tq = —ip + tt/2N, h = 1 in representation (44) for U = Un^, 
we obtain the corresponding Cartesian coordinates (xo,yo) / 0. Then 
U(x ,y ) = and applying the Euler theorem on homogeneous functions 
yields 

dU 

(51) x ■ U' x (x , yo) + yo • "^"( x o, Vo) = k*U(x 0l y ) = 0. 

Hence, the gradient VU(xq, yo) is orthogonal to the radius vector Vr(xo, yo) 
of the point (xq, yo). On the other hand, one can readily see from (44) that 

x = (//* - 1) cos - $ 

yo = (^*-l)sin(^-V 
hence we have for the polar angle of the point (xQ,yo) 

(52) o = JL_^/, modvr. 

By (46), the gradients of un and U are mutually orthogonal. Then from 
(51) we infer that vectors Vujv(:ro, yo) ano ^ Vr(xo,yo) & r e collinear. Since 
Vu N = f N {6)r k - l Vr + r k f' N (9)V8, 

and 

(V0; Vr) = 0, 

we conclude that f' N (6o) = 0. Hence, by (25) there is n € Z such that 

Tn _im 
9o = ^ = W 

which by (52) yields Nip = ^ mod it, and the theorem follows. □ □ 
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5. Algebraic iV-solutions 

In this section we settle the following question: For which rational numbers 
7 G Q such that \ j\ > 1, equation (1) does admit nontrivial (i.e. N > 2) al- 
gebraic solutions? Clearly, algebraicity of a general iV-solution is equivalent 
to that property for basic iV-solutions. Therefore, in what follows we cofine 
ourselves by only basic iV-solutions. 

Lemma 5.1. Let 7 € Q, |7| > 1 and N G N. Then inclusions £(7, N) G Q, 
A (7, N) G Q and /i(7, N) G Q are pairwise equivalent. 

Proof. It immediately follows from (27) that inclusions k G Q and A G Q 
are equivalent. By (38), (i £ Q implies both k G Q and A G Q. On the 
other hand, if k G Q then A G Q, so by (37) fi £ Q and the lemma is 
proved. □ □ 

For the further convenience, we put 7 G A, if there exists iV > 2, iV G N, 
such that un(x,v) is an algebraic function. In this case we also define 

Af("f) = {N G N : un is an algebraic function}. 

Lemma 5.2. Let 7 G A, M > 1 and N > 2. Then N G ^(7) #f the 
corresponding exponent £(7, iV) G Q. 

Proof. Let TV G ^(7). Then it follows from (44) and \fi\ > 2N - 1 > 3 that 

(x 2 + y 2 ) = h 4N - 2 [fi 2 + 1 + 2/icos(2iVT + 20)] ~ /j 47V - 2 , ^ 00 
while 

u n (x, y) = ^ i ( 2JV - 1 ) fe C os(2Afr + 0), 
where k = k(~f,N). Thus, the growth exponent of u is equal to k and it 
follows that k G Q. Now, let k(j,N) G Q. Then (45) gives a rational 
parametrization of ujf, where d is the denominator of k. Hence, un(x,u) is 
an algebraic function. □ □ 

Corollary 5.3. For 7 = 1 all the N -solutions are algebraic functions, i.e. 

M(i) = n. 

The next assertion is an easy corollary of (50) 
Lemma 5.4. Let 7 G Q. Then k^ G Q 44> k* N G Q. In particular, 

M{ 1 )=N{- 1 ). 

By Lemma 5.4, we can assume without loss of generality that 7 > 1. 
In what follows, we suppose that p and q have no common divisors. Then 
by Lemma 5.2, 7 G A if and only if there is an integer N > 2 such that 
k(j,N) G Q. By virtue of Lemma 5.1, this is equivalent to existence of a 
rational solution fi = jj > 1 of (39): 

a(A 2 - B 2 ) 



ALGEBRAIC STRUCTURE OF QUASIRADIAL SOLUTIONS 



15 



Thus, we arrive at the following diophantine equation 

(53) A 2 q - 2ABpN + q(2N - l)B 2 = 0. 

Theorem 5.5. The following assertions are equivalent 

(i) ~f = p /q£A with N £ A/"(7), N > 2; 

(ii) equation (53) has an integer solution (A, B): A, B € Z, and A > B > 

i; 

(Hi) the discriminant 

(54) D = N 2 p 2 - q 2 (2N - 1) 

is an squared integer. 

Moreover, if j £ A, 7 7^ 1, then the set ^(7) is finite and the following 
upper bounds holds 

q 2 {p 2 + 2-q 2 ) 



(55) N < ^ 

In particular, 



(56) q > 3. 
and 

(57) <? + 1 <p < - 1; 

Proof. Clearly, we have only to establish the equivalence (ii) and (iii). In 
turn, the only nontrivial implication is (iii) => (ii). 

Let (iii) be true. Then p = qj > q and D = d 2 for some d G N. Since 

D = N 2 {p 2 - q 2 ) + q 2 (N - l) 2 > 0, 

we have d > 0. Let V = A/B and consider the following associated with 
(53) quadratic equation 

(58) F(V) := V 2 q - 2pNV + q{2N - 1) = 0. 

Then (58) has two distinct rational solutions v\ and V2, vi < v-i- Since 

F(2N - 1) = -2(N - l)(2N - l)(p - q) < 0, 

we have t>2 > 2A^ — 1. Let i>2 = A/S where A and -B have no common 
divisors, and B > 0. Hence, ^4 > (2N - l)B > S and (A, B) is a desired 
solution of (ii). 

In order to establish the finiteness of J\f(j) , we make use (iii) . We have 

,50, D .(^_^) 2 + rt^!l. 

Let D = d 2 , with d € N. Then (59) yields 

q 2 

(60) d>Np-—, 

p 



16 



VLADIMIR G. TKACHEV 



while the Bernoulli inequality and (59) imply the upper following bound 
(61) d< ( Np -l)+ 



p J 2p(Np 2 - q 2 ) ' 
Since p and q have no common divisors, we can write 

Q 2 m 

— = M H 

p p 

where M > is an integer, and m € {1, 2, . . . ,p— 1} is the non-zero remain- 
der of the latter ratio. On the other hand, since m = q 2 — Mp, it follows 
that m and p have no common divisors. Thus, by virtue of Np € N, and 
strict inequalities (60) and (61) we obtain 

?V -q 2 ) i 



> -, 



2p{Np 2 -q 2 ) ' p 

which easily implies (55). 

To verify (56) we notice that the cases q = 1 and q = 2 together with (55) 
easily yield the contradiction: N < 2. 

The left inequality in (57) immediately follows from p > q. Finally, to 
prove the right hand side inequality we put as above D = d 2 , d > 1. Then 
D < N 2 p 2 , whence d < Np. Taking into account integrity of d we obtain 
d < Np — 1, or what is the same 

< {Np - l) 2 - D = 2N(q 2 -p) + l-q 2 , 

hence 

(62) q 2 -p > (q 2 - l)/2N > 0. 

Therefore, p < q 2 and by virtue of integrity of q and p we arrive at a stronger 
inequality 

q + 1 < p < q 2 ~ 1 

and the theorem is proved. □ □ 

Corollary 5.6. If 7, \j\ > 1, is an integer number then (I) has no (non- 
trivial) algebraic N -solutions. 

Proposition 5.7. The following representation holds 
f 2/V — 1 + s 2 

(63) ^1= — ^ : seQn(0,l), NGN 

I 2sN 

Proof. We notice that the discriminant D in (54) will be the a squared 
integer if and only if a positive integer valued solution (x, y) to the following 
Pell type equation 

(64) N 2 x 2 - (27V - l)y 2 = 1, 

does exist such that x > y. Here x = p/d, y = q/d and 7 is uniquely defined 
by 7 = x /y. 



ALGEBRAIC STRUCTURE OF QUASIRADIAL SOLUTIONS 17 

Then (64) can be resolved by the standard rationalization technique. Re- 
ally, let x = (l + sy)/N, whence s should be a rational number. Substituting 
the last expression in (64) gives 



2s 2N - 1 + s 



2 



V = 7T77 ; 9' x 



2N-1-S 21 N (2N - 1 - s 2 ) ' 

Thus, we have 

_2N -1 + s 2 
7 ~ 2s~N ' 

To define the admissible values of s > which provide 7 > 1 we notice 
that 

1 {s - 2N + l)(s - 1) 
7 — 1 = , 

whence s £ (0; 1) U (2N — l;+oo). We can reduce our representation to 
the interval (0; 1) only since the invariance property of (63) with respect to 
involution s -> (2N - l)/s. □ □ 

6. Maximal and minimal series 

Now we study the set A in dependence of the fractional decomposition 
7 = p/q. We suppose as before that p and q have no common divisors. 

Proposition 6.1 (Maximal series). For p/q G A the following sharp esti- 
mates hold 

(65) p < q 2 - 2 

with equality only for odd denominator q = 2s + 1, p = 4(s 2 + s) — 1, and 
N = s(s + 1), where s is a positive integer. Moreover, if the denominator 
q = 2s is even we have a stronger inequality 

a 2 -2 

(66) p < ^— — = 2s 2 - 1 



with equality iff 



g 2 ~2 q 2 - 4 

P= — ^ — , N = 



2 ' 4 

Proof. First prove (65). In view of (57), it suffices to exclude the case 
q = p 2 — 1. Assuming the contrary, and let the last equality hold. Then (62) 
and (55) together yield 

q 2 <2N<q 2 ( p2 + 2 - q2 \ 
p z 

which easily implies q 2 < 2 that contradicts to the lower bound (56). 

In order to analyze the equality case, assume that p = q 2 — 2. It follows 
then from (62) that 4iV > q 2 — 1 = p+1 and (54) can be rewritten as follows 

d 2 = N 2 p 2 - (2N - l)(p + 2) = (Np - l) 2 - (AN - 1 - p) < (Np - l) 2 
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where d > is an integer. Hence, d < (Np — 1). 
On the other hand, 

d 2 = (Np - 2) 2 + (p - 2)(2N - 1) > (Np - 2) 2 

which together with the preceding inequality implies d = Np — 1, and con- 
sequently 

p = q 2 - 2 = AN - 1. 

In particular, q must be an odd number. One can readily find that in this 
case N is an integer. Thus, the first case of the corollary is proved. 

To prove the second statement we suppose q = 2s. We have from (56) 
s > 2 and due to irreducibility of p/q we notice that p should be an odd 
number. 

By Theorem 5.5 we have for the discriminant 

D = d 2 = N 2 p 2 - As 2 (2N - 1). 

The last identity shows that d has the same parity as Np does. Thus, 
d < pN — 2. On the other hand, 

{Np - 2) 2 > d 2 = (Np - 2) 2 + A(Np - 1 - s 2 (2N - 1)) 

and therefore, 

2 s 2 - 1 
The last inequality and s > 2 yield 

n 2 

p<2s 2 = ^L. 
1 2 

Now, the inequality (66) follows from the evenness of q. 

To analyze the equality case in (66) we observe that in our notation 
d > Np — 3, otherwise we would have 

(Np - 3) 2 < d 2 = N 2 p 2 - 2(2N -l)(p+ 1), 

which implies 2N(p—2) < 7—2p, and contradiction follows. Thus, d = Np—2 
and a straightforward computation shows 

N = P -^ = s 2 -l 
2 

which completes the proof. □ □ 

Proposition 6.2 (Minimal series). For p/q G A the following lower bound 
holds 

(67) p>q + 2 

with equality iff q = 2s + 1 is an odd number, p = q + 2 and N = (q — l)/2. 
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Proof. Form (57) we have p > q + 1. Again, we argue by contradictory 
and assume that p = q + 1. Then we have the following relation for the 
discriminant 

(68) D = N 2 (q+lf -q 2 (2N -1) = d 2 , 

and d is a positive integer. Hence 

d 2 = q\N - I) 2 + 2qN 2 + TV 2 = (q(N - 1) + J^) ' - ^^jf • 



In particular, 



Af2 

d < <?(TV - 1) + — — = g (JV - 1) + AT + 1 



TV — 1 ^ v ' TV — 1 

which implies 

d< 1) +JV + 1 = di + 1. 

On the other hand 

= ( 9 (jv - 1) + TV) 2 = q 2 (N - I) 2 + 2qN(N - 1) + N 2 < d 2 , 

and it follows that d = Then by virtue of (68) we arrive at 2q = 2N—1. 

Thus, the contradictions shows that p > q + 2 and (67) is proved. 

Now assume that the equality p = q + 2 holds. Arguing as above we 
obtain the following inequality 

(69) d < (N - l)q + 27V + 2 + - < (N - l)q + 2N + 3. 
On the other hand, 

(70) d 2 = q 2 (N - l) 2 + AqN 2 + AN 2 > ((/V - l)q + 27V) 2 , 
which implies by strong inequality in (69) 

(TV - l)q + 27V + 1 < d < {N - l)q + 27V + 2. 

But it follows from (70) that d and (TV — l)q have the same parity. This 
gives exactly one choice in the last inequality 

d = (TV - l)q + 2TV + 2 

which after comparison with the definition of d in (68) implies q = 2TV + 1, 
and the required relation follows. □ □ 

The 'algebraic' constants 7 € A for small denominators q < 30 are dis- 
played in Figure 1. 
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Figure 1. q-p Diagram 



7. Concluding remarks 

We wish point out that our formula (5) can also be deduced by using the 
general representation of p-harmonic functions in the plane near their sin- 
gular points. The mentioned representation is established in [9] and [12] by- 
using the hodograph method. Our approach, nevertheless, is more direct 
and use no the quasiregularity of the complex gradient of the ^-solutions. 

The author wish to thank the anonymous reviewer for many constructive 
advices and bringing our attention to the papers [9] and [12]. 
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